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ABSTRACT 

We present early results from the JCMT Plane Survey (JPS), which is surveying the 
northern inner Galactic Plane between longitudes i = 7° and i = 63° in the 850- 
pm continuum with SCUBA-2, as part of the James Clerk Maxwell Telescope Legacy 
Survey programme. Data from the £ = 30° survey region, which contains the massive 
star-forming regions W 43 and G29.96, are analysed after approximately 40% of the ob¬ 
servations had been completed. The pixel-to-pixel noise is found to be 19mJybeam“^ 
after a smooth over the beam area, and the projected equivalent noise levels in the 
final survey are expected to be around lOmJybeam”^. An initial extraction of com¬ 
pact sources was performed using the FellWalker method resulting in the detection 
of 1029 sources above a 5-cr surface-brightness threshold. The completeness limits in 
these data are estimated to be around 0.2Jybeam“^ (peak fiux density) and 0.8 Jy 
(integrated fiux density) and are therefore probably already dominated by source con¬ 
fusion in this relatively crowded section of the survey. The fiux densities of extracted 
compact sources are consistent with those of matching detections in the ATLASGAL 
survey. We analyse the virial and evolutionary state of the detected clumps in the 
W43 star-forming complex and find that they appear younger than the Galactic-plane 
average. 
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1 INTRODUCTION 

In order to obtain a predictive understanding of star forma¬ 
tion that will be useful to all areas of astrophysics, we must 
identify and quantify the principal mechanisms responsible 
for regulating the star-format ion rate (SFR) and efficiency 
(SFE) on different spatial scales, and, potentially, those that 
determine the stellar IMF. We need to separate the effects of 
Galaxy-scale mechanisms, in particular the spiral-arm den¬ 
sity waves, from local triggering by feedback due to radia¬ 
tion, stellar winds and Hii regions. 

Until recently, the study of star formation (and espe¬ 
cially massive star formation) has been limited by biased and 
incomplete samples. For example, what controls the SFE, 
whether the IME is variable, or even which regions are typi¬ 
cal and which are extreme are all unknowns. However, these 
questions are now solvable by using large surveys to measure 
the star-format ion efficiency (SEE) and luminosity function 
(LE) of YSOs. Using Red MSX Source (RMS) survey data, 
( [Lumsden et al.|2013 ), Moore et al.| ( |2012| ) have found that 
around 70% of the increase in star-format ion rate density 
associated with spiral arms in the inner Galaxy is due to 
source crowding, rather than any physical effect of the arms 
on the star-format ion process. The remainder could be due 
to increases in SEE or in the mean luminosity of massive 
YSOs (i.e. a flattening of the MYSO luminosity function). 
However it might also be the result of residual bias in the 
source samples due to patchy star formation within spiral 
arms ( [Urquhart et al.|2Q14a ). 


Moore et al. 


( 2012| also show evidence that the mean 


molecular cloud mass decreases with Galactocentric radius, 
implying a significant change in the star-forming environ¬ 
ment with Galactic radius. The RMS survey is also begin¬ 
ning to unravel the evolutionary timescales of massive young 
stellar objects (MYSOs). Mottram et al. (2011) find that 
the accretion phase of MYSOs has a duration ranging from 
4 X 10^ years for objects with luminosities of 10^ L© to 
~ 7 X 10^ years at for those with L — 10^ L©. At lumi¬ 
nosities above 10^ L© the accretion timescales for MYSOs 
become comparable to the Kelvin-Helmholtz time. |Davies| 
et al. (2011) suggest that these timescales are most consis¬ 


tent with accretion rates predicted by turbulent core and 
competitive accretion models (e.g. Bonnell k, Bate||2006 ). 


Similarly, there has been recent progress on constrain¬ 
ing the role of spiral arms in affecting the efficiency of form¬ 
ing dense, star-forming clumps out of the larger molecu¬ 
lar clouds (the clump formation efficiency, or GEE). Using 


data from the Galactic Ring Survey (GRS: Jackson et al. 
2006) and the Bolocam Galactic Plane Survey ([Aguirre ^ 


al.|2011 ), in the range t — 37.8° to 42.5°, Eden et al. (2013)^ 


found no difference in the GEE between the arm and inter¬ 
arm regions. This implies that the increased star formation 
rate in Galactic spiral arms arises simply from the presence 
of greater amounts of gas in the arms, and is not due to an 
increase in the efficiency of forming dense clumps. 

The key data required in such investigations are pro¬ 
vided by large-scale surveys of the Galactic Plane in tracers 
of the cool, dense, molecular gas which forms the initial con¬ 
ditions for star formation. Of particular importance are sur¬ 
veys in rotational emission lines from the isotopologues of 
GO and continuum emission in the far-infrared and sub- mil¬ 
limetre bands from warm and cold dust mixed with the gas. 


The thermal sub-millimetre continuum emission from cool 
and cold dust is associated with a variety of early phases 
of star formation and is reliably optically thin, giving di¬ 
rect measurements of column density if it is assumed that 
gas and dust are well mixed. Large-area surveys at these 
wavelengths therefore provide a census of current and incip¬ 
ient star-format ion activity. Unbiased, blind, surveys are re¬ 
quired to demonstrate both where star formation is present 
and where it is absent. 

A number of major continuum surveys of the Galactic 
Plane (GP) have been completed or are in progress, which 
will allow studies such as those above to be made in an un¬ 
biased fashion, with more detail, and over Galactic scales. 
These surveys hold out the promise of providing data that 
will facilitate major steps forward in the study of star forma¬ 
tion by providing SED and luminosity information for com¬ 
plete samples of star-forming regions and individual young 
stellar objects, large enough to bin by luminosity, evolu¬ 
tionary stage and environment. Among these are: Hi-GAL 
(Molinari et al. 2010a|b ) - an open-time Key Program of 
the Herschel Space Observatory, which carried out a 5-band 
photometric imaging survey at 70, 160, 250, 350, and 500 
pm of a \b\ ^ 1° wide strip of the entire Galactic plane and 


following the Galactic warp and ATLASGAL (Schuller et 


al. 2009), the first systematic survey of the inner Galactic 
plane at 870 pm, covering 280° < £ < 60° and \b\ < 1.5° 
over most of its range. This survey was carried out using 
the LABOGA instrument at the APEX telescope |Siringo~^ 
al.|2009 GORNISH ( Hoare et al.|2012 - which has used the 


VLA to map the northern Galactic plane at 5 GHz in the 
range £ = 10° — 60° and at latitudes of \b\ < 1°. GLIMPSE 
( Benjamin et al.||2003 Ghurchwell et al.]|2009 ) - which sur¬ 
veyed a range of Galactic latitude and longitude at wave¬ 
lengths of 3.6, 4.5, 5.8, and 8.0 pm. Specifically £ — 10° —65° 
at \b\ < 1°, the inner 20° of the Galactic plane with b < 2° 
in the inner ^ = 0° — 5°, and then nine selected strips ex¬ 
tended in latitude to ±3° (±4° if they were within 2° of the 
Galactic Gentre). 

In this paper, we describe early results from the James 
Glerk Maxwell Telescope (JGMT) Galactic Plane Sur¬ 
vey (JPS), part of the JGMT Legacy Survey programme 
( Ghrysostomou|[2009 ). The results presented are from a re¬ 
gion of the survey near Galactic longitudes of ^ = 30°, which 
is well enough studied so that consistency checks can be 
made and where the source density is high and the detection 
sensitivity is near to the confusion limit and unlikely to be 
much improved upon in the final product. JPS was allocated 
450 hours in JGMT weather band 3 (0.08 ^ T 225 GHz ^ 0.12) 
to survey the region 7° ^ £ ^ 63°, \ b\ ^ 1° in the 850- 
pm continuum. Observations began in June 2012 and were 
completed in late 2014. The JPS Gonsortium consists of as¬ 
tronomers affiliated at time of joining to one of the three pre- 
2013 JGMT partner countries: UK, Ganada & the Nether¬ 
lands. JPS is part of a coordinated plan for mapping the 
Galactic Plane in the sub-mm continuum with JGMT, which 
involves a second JGMT/SGUBA-2 Legacy-survey project, 
the SGUBA-2 Ambitious Sky Survey (SASSy; Thompson et 
ah, in preparation). SASSy will observe the outer GP region 
120° < £ < 240°; \ b\ ^ 2.5° at 850 pm with sensitivity 30- 
40 mJybeam”^ rms. The SASSy Perseus extension will also 
map longitudes 60° < £ < 120°, completing the SGUBA-2 
coverage of the GP visible from JGMT to an approximately 
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constant average mass sensitivity. JPS covers the same area 
as the ATLASGAL survey ( Schuller et al.||2009 ) and Bolo- 
cam Galactic Plane Survey (BGPS; [Aguirre et al.|2011 ) but 
will be several times deeper. 

The goal of the JPS is to produce a census of, princi¬ 
pally, massive star formation across a large fraction of the 
Galaxy, complete to a mass limit of order 100 Mq (depend¬ 
ing on assumed temperature and emissivity values) at a dis¬ 
tance of ~20 kpc. Aside from the legacy value of the data, 
the immediate science aims are to investigate the earliest 
stages of the evolutionary sequence of massive protostars 
and pre-stellar objects, the variation of star-format ion effi¬ 
ciency with location in the Galaxy and the effect of envi¬ 
ronment and feedback on this and the mass distribution of 
dense, star-forming clumps. A public release of JPS data will 
accompany a future paper describing the complete survey. 


2 OBSERVING STRATEGY 


JPS uses the wide-field submillimetre-band bolometer cam¬ 
era SGUBA-2 (the Sub-mm Gommon-User Bolometer Array 
2: Holland et al.|2()T3 ) in the 850-pm band at a spatial res¬ 
olution of 14.5 arcseconds. 

JPS will sample the region of the inner GP within the 
longitude range 7° < £ < 63° and latitude \b\ ^ 0.8°, an 
area coincident with that of the northern GLIMPSE survey 


(Benjamin et al.||2003 

also MIPSGAL: 

[Garey et al.j 

12009), 

and GORNISH north ( 

Hoare et al.|2012 

jPurcell et al.| 

|2013). 


The strategy is to map six large, regularly spaced fields 
within this area, centred at intervals of 10° from £ = 10° 
(Figure]^ with a target rms sensitivity of ~10mJybeam“^ 
at 850 pm. Each of the fields covers approximately 5° in lon¬ 
gitude and 1.7° in latitude. This strategy was developed in 
response to the re-scoping of the entire JGMT Legacy Sur¬ 
vey project that took place in 2011. The choice of regularly 
spaced but spatially limited target fields allows a significant 
increase in depth over existing surveys while preserving the 
goal of producing a relatively unbiased survey and still tak¬ 
ing in several significant features of Galactic structure, par¬ 
ticularly the tangent and origin of the Scutum spiral arm at 
£ = 30° and the Sagittarius-arm tangent at I = 50°. Figure 
shows the planned extent of JPS in the context of several 
other key sub-mm surveys. 

The observing modes developed for SGUBA-2 include 
large-area mapping modes known as rotating pong patterns 
in which the telescope tracks across a square area, filling 
it in by notionally bouncing off the boundary of this area, 
in the manner of the eponymous vintage video game. Once 
a pattern is complete the map is rotated and the pattern 
repeated at the new angle. The combination of telescope 
velocity, the spacing between successive rows of the basic 
pattern and the number of rotations have been optimised 
to give the most uniform exposure-time coverage across the 


defined field. JPS uses the pong3600 mode (Bintley et al. 


2014), with 8 rotations and a telescope scan speed of 600 


arcseconds per second to produce a one-degree circular area 
of consistent exposure time, to observe each tile in the survey 
(the actual mapped area is larger than this). The total time 
required for each observation is generally between 40 and 45 
minutes. 

The survey target areas are sampled using a regular grid 


pattern of individual pong3600 tiles with minimised over¬ 
laps. Figure shows an exposure-time map of part of such 
a grid. The overlaps, in which the exposure time is roughly 
doubled, constitute about 9% of the grid area. Individual 
pong3600s have a pixel-to-pixel rms noise of 50 to 80 mJy 
per beam when reduced with 4-arcsecond pixels and each tile 
of the final survey will contain approximately seven repeats 
of these. As far as possible, even coverage of the whole tar¬ 
get area was maintained as the survey progressed, in order 
to produce relatively uniform sensitivity at any stage. 

The 850-pm survey data presented in this paper cover 
the £ = 30° field of the JPS and were observed between 
June 2012 and October 2013. The 11 tiles making up the 
field were observed on average three times each. A strategy 
of minimum elevation limits for given atmospheric opacity 
bands within the allocated range was adopted, in order to 
minimise variations in the resulting noise in each repeated 
tile. 


Flux calibrations, pointing and focus checks were made 
routinely during each observing session with observations 


of isolated observatory standard point sources (Dempsey et 
^|2013a ). The flux calibrations were typically found to be 


repeatable to within 5%. The pointing accuracy is typically 
less than 5 arcseconds rms. 

SGUBA-2 observes the 850-pm and 450-pm bands si¬ 
multaneously but the assigned weather band for JPS obser¬ 
vations ensures high sky opacity (r ~ 1.7 to 2.7 in weather 
band 3) at 450 pm. The brightest sources are detected at 
low sensitivity at 450 pm but the fluxes are not photometri¬ 
cally reliable and therefore are not currently included in the 
survey science data. 


3 DATA REDUCTION 


The data presented in this paper were reduced with the Dy¬ 
namic Iterative Map-Maker ( Ghapin et al.||2013 ), available 
as part of the Starlink SMURF package ( Jenness et al.|20TT ). 

The raw data are first flat-fielded, converting raw DAG 
units into pW and then sampled to the requested pixel size, 
in this case 4 arcseconds, cleaned of spikes and steps in the 
data, and finally bad data are flagged or weighted down. 

The iterative part of the reduction process then begins 
with the signal common to all bolometers being estimated 
and removed. An atmospheric correction is applied and the 
data are filtered to remove any remaining features on scales 
larger than the array footprint of 480 arcseconds, since these 
are likely to be artefacts. The remaining signal is then binned 
into a map of the sky using 4-arcsecond pixels. At this point 
the map made by the previous iteration (if any) is added 
back in, and the total map is sampled at the position of 
every bolometer value to create an estimate of the astro¬ 
nomical component of the signal, with any residual signal 
attributed to noise. These noise values are used to weight 
the bolometers when making maps on subsequent iterations, 
with each bolometer having multiple variance estimates cal¬ 
culated over a 15-second period. This astronomical signal is 
then removed from the original cleaned data, and the pro¬ 
cess is repeated iteratively starting from the modified data, 
until the change between maps produced by consecutive it¬ 
erations is on average less than 1% of the noise level. 

In order to prevent the growth of low-frequency arte- 
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Figure 1. JPS projected area coverage on the sky (green contours) overlaid on the Planck dust opacity map (colour scale; |Planck| 
|Collaboration|2014[ | . The apparent size difference between the six JPS patches is due to the tangent projection of the Galactic coordinate 
system. The zero of Declination passes through i ~ 32.9°. The axis scale is degrees. 


facts, an extra constraint is added by forcing all background 
pixels in the map to be zero. Background pixels are taken 
to be those with a signal-to-noise ratio less than 3.0, plus 
any other contiguously attached pixels down to a signal-to- 
noise ratio of 2.0. No such constraint is applied on the very 
last iteration, thus allowing non-zero values to appear in 
the background regions of the final map. This means that 
the background values within the final map are the result 
of a single iteration, while the astronomical signal has been 
arrived at by many iterations. 

The present data have been reduced as individual pong 
tiles onto a 4-arcsecond pixel grid and then mosaiced to con¬ 



Figure 2. The area of the Galaxy covered by JPS (green seg¬ 
ments). Also shown, for comparison, is the longitude coverage of 
ATLASGAL (yellow), BGPS (blue) and SASSy (pink). The back¬ 
ground image is the sketch of the Milky Way by Robert Hurt of 
the Spitzer Science Genter, made in collaboration with Robert 
Benjamin. 


struct the whole I — 30° field map. Prior to creating the mo¬ 
saic, noisy tile edges were removed by cropping each observa¬ 
tion to a 40 arcminute radius. The combined map was then 
lightly smoothed using a Gaussian kernel with a = 1.8 pixels 
(7.2 arcsec), approximating the half-width at half maximum 
of the telescope beam at 850 [im. The standard observatory- 
determined calibration factor of 537 Jy beam“^ pW“^ was 
applied to the data. 
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4 RESULTS 


Figure shows all of the reduced i = 30° field data. Well- 
known star-forming regions W 43 and G 29.96 can be iden- 
tihed at positions i = 30.815, b = —0.057 and i = 29.954, 
b = —0.018, respectively. Figure shows a close-up of the 
same data in the 1-degree region around W 43. 

The distribution of all pixel values in the I = 30° held 
is shown in Figure along with a Gaussian ht to the distri¬ 
bution (—0.1 to 0.1 Jybeam“^). The rms noise in the data 
is estimated from the latter to be arms = 19mJybeam“^. 
It can be seen from the negative pixel values that there is 
a slight excess noise in the wings of the distribution above 
that of the assumed Gaussian. The larger excess on the pos¬ 
itive side is due to signal from real sources. Since the data 
set used here consists of around 40 per cent of the eventual 
JPS data, the projected rms sensitivity in the complete sur¬ 
vey, in regions where confusion is not dominant, is therefore 
likely to be around 12mJybeam“^. This can be compared 
to the ATLASGAL rms sensitivity of 50 - 80 mJy. 


60 

50 

40 

30 

20 

10 

0 

Figure 3. Part of the survey grid pattern consisting of overlap¬ 
ping 1-degree circular pong tiles. The colour scale indicates the 
spatial distribution of relative exposure time (in arbitrary units) 
assuming identical observing conditions and the same number of 
repeats for each tile. 



4.1 J=3-2 line contamination 


With a central frequency of 355 GHz and half-power width 
around 35 GHz, the continuum bandpass at 850 pm will 
suffer contamination from spectral-line e mission, predomi- 
nantly by ^^GO (J=3-2) at 345.796GHz (Johnstone, Boon- 
man van Dishoeck||200^ . To quantify the degree of con¬ 
tamination, we reprocess the SGUBA-2 data with the GO 
contribution removed during the reduction. 

GO data for this region were taken from the GO High 
Resolu tion Survey (GOHRS: [Dempsey, Thomas fc Gurrie 


2013b). GOHRS has mapped the Galactic plane in GO 


(J=3-2) in the Galactic longitude range £ = 10° to ^ = 
55° using HARP on the JGMT and its spatial resolution is 
therefore essentially identical to that of the SGUBA-2 data. 
GOHRS only extends to Galactic latitudes b = ±0.25° in the 
£ = 30° region, so only this central strip can be considered 
in this analysis. 

The technique for removing the GO contamination in¬ 
volves running the mapmaker as described in Section 3 but 
with the GO map, converted to SGUBA-2 output power 
units, supplied as a negative fake source. This method has 
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Figure 4. JPS 850-pm surface-brightness distribution in the £ = 30° field as of October 2013 with about 40 per cent of the data acquired. 
Individual circular pong tiles have been cropped at a radius of around 0.65°. The bright sources near longitudes of 30.0° and 30.8° are 
the well-known star-forming regions G29.96 and W43, respectively. The colour scale range is logarithmic from zero to 23Jybeam“^ 
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Figure 5. 850-|xm surface-brightness distribution in the ~l-degree square area around the W 43 star-forming region. The data have 
been lightly smoothed with a cr = 7.2" Gaussian kernel, equivalent to the 850-|J.m JCMT beam, slightly degrading the native spatial 
resolution. 


previously been used by Drabek et al.| ( |2012| ), [Hatchell"^ 


al. (2013) and Sadavoy et al. (2013) to investigate CO con¬ 


tamination. To convert to physical units (pW), we begin by 
multiplying the CO intensities (in Kkms“^) by the forward¬ 
scattering and spillover efficiency (0.71) and by a molecular- 
line conversion factor. This factor is calculated from the at¬ 
mospheric precipitable water vapour (P) and scales with 
opacity at 225 GHz using the following equation, described 
in full in Parsons et al. (in preparation): 


C850 = 0.574 1 7 + 0.115089 P - 0.048470 P^ + 0.010893 P^ 
— 0.000856 P^ mJybeam~^/(Kkms~^) 

The final step converts the CO data from mJybeam”^ 
into negative pW by dividing by the standard 850-qm cali¬ 
bration factor (—537000 mJybeam”^ pW“^). Once the raw 
data have been reduced with the CO contribution removed, 
the resulting map is compared to the uncorrected SCUBA-2 


data, applying a mask based on a 3-cr cut in the signal-to- 
noise map of the original reduction to both maps. 


For the bulk of the data, the levels of CO contamination 
range from zero to ~30% (see Figureand the majority of 
values are only a few per cent. This is consistent with values 
reported by Schuller et al.| ( |20Q9| who quote results from 


Wyrowski et al. (2006) and with those found by Drabek et 


al. (2012). The distribution is not random but varies co¬ 


herently between and within sources, with clearly resolved 
gradients that are sometimes radial and sometimes longitu¬ 
dinal. In general, the highest levels of contamination appear 
towards the edges of continuum sources, where the ^^CO 
optical depth is likely to be least. Detailed analysis of the 
spatial variation in CO contamination is beyond the scope 
of this paper and will be addressed in a future study; how¬ 
ever, for the brightest sources with a SCUBA-2 flux between 
1000 and 1500mJy beam“^, the contamination levels have 
a mode of 2% and do not exceed 10%. It is therefore clear 





































that the level of contamination is relatively low and unlikely 
to affect significantly the extracted dust column densities. 

Since the JPS has a wider latitude range (| ^ | ^ 0.85°) 
than COHRS, a statistical correction would be desirable 
where no COHRS data are available. However, we found 
no correlation between the level of contamination and the 
surface brightness in the continuum map and so such an 
extrapolation is unlikely to be possible. 


Another potential source of contamination of the ther¬ 
mal dust emission at sub-millimetre wavelengths is that from 
the free-free continuum. Obviously, this would be most se¬ 
vere for luminous compact Hll regions, such as those in W43. 


Schuller et al. (2009), using results from Motte et al. (2003), 
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estimate this contribution to be less than ~20 per cent in 
most cases. 


4.2 Detection of compact sources 


To identify compact sources within the data, we have used 



gorithms, FW performs best when the background noise is 



Figure 6. The distribution of all pixel values in the I = 30° 
field (black histogram) and the result of a Gaussian fit to the 
distribution (—0.1 to O.lJybeam”^; dashed red curve). The sen¬ 
sitivity of the current data is estimated from the latter fit to be 
Crms = 19mJybeam“^. 



Figure 7. Spatial distribution of the J=3-2 emission line contamination in the region around W43. The colour scale shows the 

value of (5(850) — 5(^^CO))/5(850) and the range is from 0.7 (darkest) to 1.0 (lightest). 
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uniformly distributed. We have therefore used the weight 
map produced by the data-reduction processing to estimate 
the rms noise distribution in the image and so produced a 
signal-to-noise ratio (SNR) map of the whole region. The 
noise in this SNR map is effectively constant, even towards 
the edges where the rms in the surface-brightness map in¬ 
creases significantly. 

The SNR map was used as the initial input into FW and 
all sources above a specihed SNR threshold were identified 
with the task FiNDCLUMPS. The FW algorithm assigns pixels 
to clumps by ascending towards emission peaks from each 
pixel above the threshold via the steepest gradient. It then 
creates a mask of the input map in which pixel values either 
identify the source to which the pixel has been assigned or 
are set to a negative value if unassigned. This mask is then 
used as input to the task Extractclumps that extracts 
peak and integrated flux values from the original emission 
map for each source identified. 

We set a detection threshold of 3cr and required that 
the source consisted of >7 pixels; this is the number of pix¬ 
els we would expect above the detection threshold for an 
unresolved 5-cr source. The sensitivity of JPS data to the 
extended, filamentary structures found to be ubiquitous in 
Herschel Hi-GAL images ( Molinari et al.|2010b ) has not yet 
been tested and so, for the purposes of this paper, we have 
rejected sources with aspect ratios larger than 5. 

Figurej^ displays the distribution of 850-pm surface 
brightness in a small subregion of the I = 30° field, around 
the prominent W43 star-forming region, and the sources ex¬ 
tracted from the data by FW. Most of the sources that are 
identifiable by eye have been recovered, although the de¬ 
tection of sources below 5cr is relatively unreliable, further 
justifying a threshold at this level for inclusion in our prelim¬ 
inary compact source catalogue. The latter consists of 1029 
objects and is listed in Table The set of FW parameter 
settings used is given in Table 


4.3 Flux calibration and photometric corrections 

FellWalker produces peak flux density values taken sim¬ 
ply from the highest pixel value associated with each ex¬ 
tracted source. Integrated flux densities are calculated by 
summing the values of all associated pixels in the source 
mask and dividing the result by the so-called beam inte¬ 
gral, which is the number of pixels per beam. The latter 
is determined from the ratio of integrated to peak surface 
brightness of a point source. In this case, the observations 
of Neptune, reduced in the same manner as the science data 
were used, producing a beam integral value of 17.5, which is 
significantly larger than that expected for a Gaussian beam 
(~14), although the FWHM of the Neptune image is consis¬ 
tent with the latter. This is due to significant power in the 
wings of the JGMT beam, as discussed by |Dempsey et al.| 
( 2013^ . 

In addition to this, the FW method, in common with 
similar algorithms, effectively uses variable aperture pho¬ 
tometry to produce integrated flux densities, with the aper¬ 
ture size set by the number of associated pixels above the in¬ 
clusion threshold in the source masks. This number is SNR- 
dependent and so requires the application of an aperture 
correction. This has been done by converting the area in pix¬ 
els {A) of each source to an effective source circular radius 


r = d\JAji:^ where d is the angular pixel scale, and then us¬ 
ing the aperture correction values tabulated by [Dempsey et| 


al. (2013a). The latter were checked using our own observa¬ 


tion of Neptune, and found to be consistent to within 1 per 
cent. Values intermediate to those tabulated were interpo¬ 
lated using a polynomial fit (see Figure]^. Note that these 
corrections are only approximate, since they take no account 
of departures from circular symmetry. However, they only 
signihcantly affect small sources, those with sizes only a lit¬ 
tle larger than the beam, and these are likely to be circu¬ 
lar. Most irregular sources are extended with respect to the 
beam and for sources with diameters > 60 arcseconds the 
correction is of order 10 per cent or less. It is also the case 
that very compact and unresolved sources are relatively rare. 
The error beam contribution is estimated to have an ampli¬ 
tude of 2 per cent (Dempsey et al. 2013a), which is much 
less than the other uncertainties. 


As a check on the photometry of the sources extracted 
by FellWalker, we have compared the peak flux densi¬ 
ties of the JPS compact sources to those of spatially corre¬ 


sponding detections in the ATLASGAL survey (Gontreras 
et al.|[20T3 ). We searched the ATLASGAL compact -source 


catalogue ( Urquhart et al.|2014b ) for the closest spatial co¬ 
incidence within 10 arcseconds of each JPS detection and 
found 484 matching sources. Figure 10 (upper panel) shows 
the comparison of peak flux densities for this initial matched 
list. It can be seen that the JPS peak flux densities are sys¬ 
tematically lower than those of ATLASGAL, with a scat¬ 
ter that increases towards lower values. Much of this sys¬ 
tematic difference can be attributed to the fact that most 
sources in the ATLASGAL catalogue were found to be at 
least somewhat extended. Given the larger APEX beam (19 
arcsec, compared to around 14.5 arcsec for JGMT) at these 
wavelengths, we expect peak flux densities to be larger in 
the former data. To correct for this, we repeated the JPS 
source extraction after convolving the data to the ATLAS¬ 
GAL spatial resolution and scaling the signal level so as to 
restore the (equality of) peak and integrated flux densities 
of point sources. The lower panel of Eigure shows the re¬ 
sult, this time with 500 positional matches and a significant 
improvement in the flux matching. The mean ratio of peak 
flux densities in the two catalogues for these latter matched 
sources is aSsso(JPS)/aSsto(ATLAS) = 0.93±0.01 and a stan¬ 
dard deviation of 0.12. Despite being observed with different 
telescopes and detectors and using unrelated techniques, re¬ 
duced using different procedures and the source extraction 
being done using independent methods (ATLASGAL uses 
Sextractor), the fluxes in the two catalogues are in very 
good agreement. This is, therefore, an important consistency 
check on the reliability of the results of both catalogues. 


The scatter in the correlation widens at lower flux den¬ 
sities as might be expected where the signal-to-noise ratio in 
the ATLASGAL data becomes low. However, there is also 
a deviation from the correlation at the low end, in the form 
of a tail with high values of aSsso(JPS)/aSsto(ATLAS). This 
tail may be the result of a bias related to the flux-boosting 
bias found in the JPS data (see Section 4.5). Without it, 
the correspondence between the two flux scales becomes al¬ 
most perfect. Taking only JPS fluxes above 0.5 Jybeam”^, 
the mean peak-flux-density ratio in the matched sample is 
1.001 ± 0.001 and the median value is 1.012. 
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Figure 8. Example of the source-extraction results in a small portion of the i = 30° region. The greyscale image shows the 850-pm 
emission and the red and blue ellipses show the approximate sizes and orientations of the sources above and below 5(7. The grey contours 
start at Scr and increase in steps of 0.17, 0.8, 2.75, 6.81, 13.85 Jybeam“^. The green filled circle in the lower-left corner indicates the 
FWHM size of the JCMT beam at 850 pm. 



Q. 

< 


0.2 


0.01_,_,_,_I_,_,_,_I_,_,_,_I_I_I_^_I_^^_I_I_,_,_!_ 

0 20 40 60 80 100 120 

Aperture Diameter (arcsec) 

Figure 9. Aperture corrections determined by [Dempsey et al.| 
( |2013a| (filled circles) and re-normalised to a value of 100% at 
the largest aperture, then fitted with a fifth-order polynomial fit 
(dashed line), from which corrections were applied to the inte¬ 
grated fluxes produced by Fell Walker. 


4.4 Completeness tests 

To test the efficiency of the source-extraction algorithm as a 
function of source flux we injected artificial compact sources 
onto a field of uniform noise, the latter distributed as a Gaus¬ 


sian with a standard deviation of ~20 mJy beam“^, similar 
to that in the data. The source profiles were Gaussian and 
circularly symmetrical with a FWHM of 5 pixels (20 arc- 
sec). FellWalker was then used with the same settings 
as in the survey data to recover the artificial sources. The 
rate of detections as a function of input flux density then 
allows us to benchmark the performance of the extraction 
algorithm in this idealised case. 

The CUPID task Makeclumps was used to generate a 
catalogue of artificial objects with a uniform distribution in 
i, b. The peak fluxes of the injected objects were uniformly 
distributed between 0.02 and 0.5 Jybeam”^ and the source 
profiles were Gaussian with FWHM of 5 pixels (20''). We 
used the size of the survey field as a template and injected a 
similar number of sources as was detected by FellWalker 
in the data. FigurepT] shows the number of sources recov¬ 
ered as a function of input peak flux and that the source 
extraction algorithm is approximately 95% complete above 
80mJybeam“^ (^dcr). 

F igurep^ compares the values of the original and recov¬ 
ered peak flux densities of the injected sources. These are 
very well correlated, showing that, in general, source fluxes 
are recovered accurately. However, the recovered fluxes for 
the weaker sources are systematically higher than their in¬ 
jected values. The reason for this is that FellWalker sets 
the peak flux to be the value of the brightest pixel in the 
source and, since the injected sources are all slightly ex¬ 
tended, there is often a pixel with a positive noise deviation 
near the peak that adds constructively to the measured peak 
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Table 1. The JPS i = 30° preliminary compact source catalogue. The columns are as follows: (1) name derived from Galactic 
coordinates of the maximum intensity in the source; (2)-(3) Galactic coordinates of maximum intensity in the catalogue source; (4)-(5) 
Galactic coordinates of emission centroid; (6)-(8) semi-major and semi-minor size and source position angle measured anti-clockwise 
from Galactic north; (9) effective radius of source, i.e., i^eff = a/(A/tt), where A is the area of the source above the threshold; (10)-(13) 
peak and integrated flux densities and their associated uncertainties; (14) signal to noise ratio (SNR). 


Name 

(1) 

^max 

n 

(2) 

^max 

n 

(3) 

i 

(°) 

(4) 

b 

n 

(5) 

e’maj 

(") 

(6) 

e’min 

(") 

(7) 

PA 

(°) 

(8) 

T^eff 

(") 

(9) 

'Speak A*Spgak 

(Jy beam“^) 
(10) (11) 

*8'int '^'8'int 

(Jy) 

(12) (13) 

SNR 

(14) 

G026.830-00.208 

26.830 

-0.208 

26.830 

-0.206 

7 

6 

259 

16 

0.38 

0.08 

0.64 

0.13 

5.6 

G026.865-00.276 

26.865 

-0.276 

26.867 

-0.276 

12 

11 

112 

26 

0.36 

0.06 

1.00 

0.17 

6.8 

G026.957-00.077 

26.957 

-0.077 

26.956 

-0.076 

10 

8 

224 

26 

0.77 

0.08 

1.64 

0.17 

15.2 

G026.960-00.309 

26.960 

-0.309 

26.962 

-0.308 

14 

7 

177 

23 

0.20 

0.04 

0.67 

0.14 

5.0 

G027.000-00.298 

27.000 

-0.298 

27.001 

-0.297 

24 

11 

167 

38 

0.69 

0.07 

1.89 

0.20 

15.1 

G027.011-00.040 

27.011 

-0.040 

27.012 

-0.039 

11 

11 

151 

30 

0.47 

0.06 

1.15 

0.14 

10.4 

G027.019-00.166 

27.019 

-0.166 

27.021 

-0.167 

13 

6 

192 

21 

0.20 

0.04 

0.51 

0.11 

5.1 

G027.037-00.171 

27.037 

-0.171 

27.035 

-0.165 

22 

12 

269 

31 

0.22 

0.04 

0.98 

0.19 

5.6 

G027.066+00.013 

27.066 

+0.013 

27.078 

+0.015 

32 

16 

163 

44 

0.25 

0.05 

1.76 

0.35 

5.4 

G027.086-00.564 

27.086 

-0.564 

27.086 

-0.563 

9 

8 

149 

22 

0.39 

0.05 

0.92 

0.13 

8.7 

G027.100-00.559 

27.100 

-0.559 

27.098 

-0.557 

17 

8 

186 

30 

0.57 

0.06 

1.48 

0.17 

12.9 

G027.173 -00.008 

27.173 

-0.008 

27.173 

-0.005 

14 

10 

210 

27 

0.27 

0.05 

0.83 

0.14 

6.5 

G027.178 -00.105 

27.178 

-0.105 

27.176 

-0.105 

12 

9 

166 

27 

0.43 

0.06 

1.29 

0.17 

9.4 

G027.186-00.082 

27.186 

-0.082 

27.186 

-0.079 

29 

16 

268 

55 

3.07 

0.25 

6.19 

0.50 

67.8 

G027.187-00.148 

27.187 

-0.148 

27.187 

-0.148 

20 

11 

213 

33 

0.23 

0.05 

1.15 

0.23 

5.6 

G027.207-00.100 

27.207 

-0.100 

27.208 

-0.098 

16 

15 

214 

35 

0.44 

0.06 

2.09 

0.28 

9.6 

G027.221+00.136 

27.221 

+0.136 

27.221 

+0.136 

16 

8 

188 

30 

0.84 

0.09 

1.70 

0.18 

14.3 

G027.248+00.108 

27.248 

+0.108 

27.245 

+0.108 

16 

14 

189 

34 

0.53 

0.07 

1.64 

0.21 

10.4 

G027.256+00.135 

27.256 

+0.135 

27.255 

+0.136 

10 

7 

130 

22 

0.68 

0.08 

1.34 

0.15 

11.9 

G027.279+00.145 

27.279 

+0.145 

27.272 

+0.147 

30 

12 

193 

46 

1.66 

0.15 

6.06 

0.53 

28.5 


Note: Only a small portion of the data is provided here, the full table is only available in electronic form at the GDS via anonymous ftp 
to cdsarc.u-strasbg.fr (130.79.125.5) or via http://cdsweb.u-strasbg.fr/cgi-bin/qcat7J/ 


flux, creating a small bias. This results in the average mea¬ 
sured flux density being approximately la larger than ex¬ 
pected for the weaker sources. This systematic bias to higher 
flux values for weaker sources is well fitted by a second or¬ 
der polynomial function (Figure p^ and the recovered real 
source flux densities can be corrected accordingly. 


The polynomial coefficients are: *S'iy(corr) = —0.0374 + 
1.1875'i.(orig)—0.2655'i.(orig)^, where 5'i.(corr) and 5'i.(orig) 
are the corrected and original peak flux values (in 
Jybeam”^), respectively. We apply this correction to all 
sources with peak flux densities below 0.3 Jy beam“^. Above 
this value the correction is not necessary as it is only a few 


Table 2. Fell Walker parameter settings used in the source 
extraction 


percent of the flux value, significantly less than the ^ 5 per¬ 
cent flux calibration uncertainty in 850-qm SCUBA-2 data. 


4.5 Flux distributions (peak and integrated) 

Figureshows the distribution of peak flux densities in the 
sample of compact sources extracted from the ^ = 30 field, 
as described above and listed in Table 1. Similarly, Figure 
shows the distribution of integrated flux densities. Both of 
these are compared, without corrections for dihering beam 
size, to the equivalent distributions in the sources extracted 


from the ATLASGAL compact-source catalogue (Contreras 


|et al.|20T3l [Urquhart et al.|2014b| note that a later ATLAS- 


FELLWALKER.ALLOWEDGE 0 

FELLWALKER.CLEANITER 5 

FELLWALKER.FLATSLOPE 1 

FELLWALKER.FWHMBEAM 1 

FELLWALKER.MAXBAD 0.05 

FELLWALKER.MAXJUMP 3 

FELLWALKER.MINDIP 1.5 

FELLWALKER.MINHEIGHT 3 

FELLWALKER.MINPIX 7 

FELLWALKER. NOISE 1 

FELLWALKER.RMS 1 

FELLWALKER. SHAPE ellipse 













The JCMT Plane Survey 11 


GAL catalogue is available: Csengeri et al.|2014 | in the same 
area of sky. Both figures include a power-law fit to the JPS 
data above the turnovers in each case. 


The distributions in Figure are very similar, with 
much the same slope that approximates well to a power law 
in both data sets. The main difference is that the ATLAS- 
GAL distribution turns over at a peak flux density about 
five times higher than that of the JPS sample and has a 
slight excess of brighter sources. In contrast, the integrated 
flux-density distributions in Figure are significantly dif¬ 
ferent. The ATLASGAL sample distribution is shifted to 
values that are around a factor of two higher than for JPS 
and appears to be less well represented by a power law. The 
latter difference is likely to be the result of source merging 
in the larger APEX telescope beam FWHM). 

The turnover in the peak flux-density distribution, 
which can normally be assumed to mark the effects of incom¬ 
pleteness in the sample, is significantly higher than the limit 
calculated from the tests on artificial data described above. 
This is likely to be because of confusion due to source crowd¬ 


ing and perhaps non-uniform noise, effects that are likely to 
be position-dependent within the survey. If this is the case, 
the true completeness limit in the i = 30° field may be 
as high as 0.1-0.2 Jybeam“^ (which may be compared to 
the correspon ding limit in the ATL ASGAL survey at 0.3- 
0.5 Jybeam“^, Gontreras et al. 2013). Estimating again from 


the turnovers in the distributions, the JPS integrated source 
flux densities appear complete above ~0.8 Jy while ATLAS¬ 
GAL limit is around 2-3 Jy. More rigorous tests will be con¬ 
ducted to quantify completeness when the entire JPS data 
set is available. 


4.6 The Mass Distribution of Clumps in W 43 

The W 43 Giant Molecular Gomplex is the most distinctive 
feature in the JPS / = 30° observations and dominates the 
emission within the region. This massive-star-forming com¬ 
plex is thought to be at the end of the Galactic Long Bar 
and molecular clouds associated with it cover longitudes of 




Figure 10. Comparison of the peak flux densities of extracted 
JPS sources to the corresponding ATLASGAL compact sources 
- top: unsmoothed data with 484 matches; bottom: with JPS 
data convolved to the resolution of ATLASGAL, producing 500 
matches. 
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I = 29-32 and latitudes of | 6 1° (Nguyen-Luong et al. 

2011) and are found at velocities of 80-110 kms“^. 


By extracting spectra at the positions of the catalogued 
JPS sources from J = 3 ^ 2 data cubes in the 

CO Heterodyne Inner Milky-Way Plane Survey (CHIMPS: 
Rigby et ah, in pr eparation), which are briefly described in 


Eden et al. 


(2012), or in CO J = 1 ^ 0 from the GRS, we 


can identify JPS sources that are in the l,b,v range associ¬ 
ated with the W 43 complex. Because of the latitude range 
of the ^^CO J = 3 ^ 2 survey, GRS data were used for 
sources with I 6 I > 0.5°. 


676 JPS sources were found in the appropriate I and 
b range for W 43, 47 of which required the extraction of 
GRS spectra to identify their velocities. 450 of these have 


velocities within the W 43 limits and so are considered to be 
associated with the complex. 

The masses of the sources were calculated using the 
standard optically thin conversion. 


M = 




KjyBj^{Td) 


( 1 ) 


Taking the mass absor ption coefficient at A = 850 pm to be 
/^850 = 0.001 m^ kg“^ ( Mitchell et al.||2001 ), and assuming 
a distance of 5.49kpc ( Zhang et al.||2014 ) and a single dust 
temperature of 18 K for all sources, the conversion becomes 
M/Mq — 225 X iSsso/Jy- 

We use the derived masses to plot the clump-mass func¬ 
tion (CMF) of the W43 complex and this is presented in 



0.00 


0.05 0.10 0.15 

Peak Flux (Jy Beam'J 


0.20 


Figure 11. The ratio of recovered to injected source counts as 
a function of input peak flux density. The dashed vertical lines 
indicate (from left to right) the Icr, 3f7, 5(7, and 7(7 noise levels. 



Figure 12. A comparison of recovered peak flux densities to the 
input values for the artificial sources used in the completeness 
tests. The solid blue line represents equality and the dashed red 
curve shows the result of a polynomial fit to the data. 
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Figure fl^ Assuming the CMF to be a power law above 
the observed turnover (log^o(^/^o) — 2.4) of the form 
AN/AM oc a least-squares ht to the CMF gives an in¬ 

dex of a = —1.87 zb 0.05. This is consistent with the results 
of |Urquhart et ah (2014c), who found values in the range 
— 1.85 to —1.91 for ATLAS GAL clumps that house trac¬ 
ers of massive-star formation, across different evolutionary 
stages and distributed over a large section of the Galactic 
Plane. 


The slope of the distribution is also consistent with 
those found for Galactic giant molecular clouds (GMCs), 
although the mass range is very different. The spectrum of 


GMC masses greater than (log^Q(M/M 0 ) = 4) can also be 
htted by a power law, wi th studies by [Solomon et al. ( 1987), 
[Heyer et al. (2001) and |Roman-Duval et al.| ( |2010 ), giving 
— 1.50 zb 0.36, —1.80 zb 0.03 and —2.02 zb 0.11, respectively 
(the Roman-Duval et al. value has been re-calculated for 
AN/AM, rather than A A/A log M, as published). 


The fact that the mass functions of the clump and cloud 
populations are indistinguishable, both across the Galactic 
Plane and in a specihc region as in this study, suggests that 
the same process, whether physical or stochastic, is deter¬ 
mining both distributions. Measurements of the fraction of 
molecular gas within clouds that has been converted into 



Figure 13. Distribution of peak flux densities of all extracted 
compact sources in the ^ = 30 held (grey shaded histogram) 
with the equivalent ATLASGAL distribution in the same pro¬ 
jected area of sky (blue outline histogram). The red line indicates 
the linear least-squares fit to the data above the completeness 
turnover. 



Figure 14. Distribution of integrated flux densities of the ex¬ 
tracted compact sources in the ^ = 30 held (grey shaded his¬ 
togram) with the equivalent ATLASGAL distribution in the same 
projected area of sky (blue outline histogram) 
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1.0 

Virial Parameter (a) 

Figure 16. Distribution of calculated virial ratios for the dense 
clumps identified with the W43 complex. 


dense clumps (the ‘clump-formation efficiency’ or CFE), us¬ 



ing the Bolocam Galactic Plane Survey (Eden et al. 


2012 


2013 Battisti Heyer]|2014 ) found a fairly constant value script ion of Euller Myers (1992) 


Figure 17. Dependence of virial ratio on clump mass in the dense 
clumps identified with the W43 complex; The linear fit shown by 
the dashed line is —0.42 ± 0.08 


velocity dispersion of the total column of gas. Using the pre- 

the average line-of-sight 


between 5 and 15 per cent when averaged over samples of 
clouds. This, together with a consistent power-law index for 
the mass functions of clouds and clumps, suggests that, out¬ 
side of the Galactic Gentre regions, the GEE is not sensitive 
to the local environment on kiloparsec scales (large varia¬ 
tions in GEE from cloud to cloud are found on smaller scales, 
Eden et al.|2Q12 ). 


however; see 


4.7 Star Formation in the W43 sample 

4-7.1 Virial State 

To calculate the virial masses of the W43 sources, the mea¬ 
sured ^^GO line width was corrected to estimate the average 


velocity dispersion is 


= Ai;eorr + 8 In 2 X 


feTkin 

ruH 


figa 


2 

fico 


77^ - b (2) 


where At'corr is the observed GO line width corrected for the 
spectrometer resolution, k is the Boltzmann constant, Tkin 
is the gas kinetic temperature, taken to be 12 K from the 
average excitation te mperatures found for W 43 clouds from 
^^GO J=2-l data by Garlhofi et al. (2013). The values in 


the latter work are derived from a single transition and are 
therefore lower limits. The choice is not critical, however, as 
in any case the thermal line-width correction to the velocity 
dispersion in Equation[^is small (no more than 10 per cent). 



10 ^ 10 ^ 10 ^ 10 ^ 10 ^ 


Clump Mass (Mg) 

Figure 15. The mass function of JPS compact sources associ¬ 
ated with clouds in the W43 complex. The red line indicates the 
least-squares linear fit to the data and the error bars indicate the 
Poisson statistics. 
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The mean masses of molecular gas and of carbon monoxide 
are |j,gas = 2.33 and |J,co = 28, respectively. The virial mass 
for each clump can then be calculated from: 


Mvir 

Me 


= 161 I -) 
pey 


AUa' 


kms“^ 


(3) 


( Urquhart et al.|2014c ) where R is the half-power radius of 
the clump, evaluated from the geometric mean of cTmaj and 
C^min from Table [T] 

The distribution of the virial ratios (a = Mvir/Mdump 
of the W43 clumps is shown in Figure The distribution 
peaks just below a = 2 with a logarithmic mean of —0.080 
and median of —0.236. 

It is a commonly observed result that virial ratios tend 
to cluster around a value within a factor of two of unity. Al¬ 
lowing for scatter due to random error, this is usually taken 
to indicate that the majority of the clumps can be consid¬ 
ered to be gravitationally bound and are likely to be either 
pre-stellar or protostellar. It should be noted, however, that 
such results are probably subject to selection effects, espe¬ 
cially that arising from the lack of stable equilibria between 
internal pressure and gravity, and the consequent short-lived 
nature of heavily gravitationally sub- or super-critical clump 
material that is observable in the 850-pm continuum. 


Kauffmann et al. (2013) (originally Bertoldi & McKee 


1992) define the critical value for the virial ratio as a = 2 


for sources which are not supported by a magnetic field and, 
as a result, sources with a < 2 are gravitationally unstable 
and require the presence of a strong magnetic field to provide 
support. This is the case for 93 per cent of the W43 sources. 

A strong relationship of decreasing a with increasing 
clump mass is found, similar to that seen in the results of 


Kauffmann et al. (2013) and by Urquhart et al. (2014c) for 


ATLASGAL clumps associated with tracers of massive star 
formation, and is shown in Figure pT A fit to the data finds a 


slope of —0.42±0.06, consistent with Urquhart et al. (2014c), 
who find —0.53 ± 0.16. This may again be the result of se¬ 
lection effects, particularly those associated with the mixing 
of gas and dust tracers in the analysis. Taken at face value, 
however, this relationship suggests that the clumps in the 
sample are not quite Larson-like. Such objects would have 
Av oc and M oc which would give a ~ const. The 
observed relation requires either a steeper dependence of M 
on R or Av to be virtually constant, or both. 


4-7.2 Evolutionary State 

W43 is a site of recent and ongoing star formation. The scale 
of the current star formation activity can be determined by 
JPS clumps that are coincident with 70-pm Hi-GAL com¬ 
pact sources, the presence of which is taken to be a reliable 
indicator of star formation (e.g. Faimali et al. 2012). 

A total of 287 70-pm sources were found to be coinci¬ 
dent with 233 clumps from the W43 complex sample, with 
the majority of JPS clumps found to be coincident with a 
single Hi-GAL source (185 clumps). The largest number of 
Hi-GAL sources found to be associated with a single JPS 
clump was 3, which happened on 5 occasions. 

The sample of massive clumps in W43 is complete above 
~100 Mq and 440 of the whole sample of 450 have masses 
above this value. If we assume that the star-format ion rate 


is constant, the relative statistical lifetimes of the pre-stellar 
and star-forming clumps above this threshold is given by the 
relative numbers in the two subsets. By taking the lifetime of 
the clumps to be 0.5 Myr (e.g. Ginsburg et al.|2012 ), we find 
that the relative lifetimes of the two phases are 0.26 x 10® yr 
and 0.24 x 10® yr for the star-forming and non-star-forming 
clumps, respectively. It is no surprise that the lifetimes of 
the two phases are very similar as the starless lifetime of an 
eventual star-forming clump has been estimated at less than 
half a Myr (e.g. Ginsburg et al.|[2012 ) and the timescale of 
the IR-bright stage of massive young stellar objects ranges 
from 7 X 10^ yr to 4 x 10® y r, for sources with lum inosities 
between 10^ L© and 10® L© ( Mottram et al.||2011 ), respec¬ 
tively. 

In a sample-based study, the luminosity of the YSOs 
detected at 70 pm can be usefully estimated from the mean 
ratio of the 70-pm flux to integrated infrared YSO luminos¬ 
ity obtained from the sample of |Eden et al. (2015). This 
ratio was measured to be log 3 ^Q(Fboi/FVo |xm) = 1.96 with a 
standard deviation of 0.56. The resulting luminosities can 
then be used to calculate the ratio L/M for the individual 
clumps. The resulting distribution oi LjM values is shown 
in Figure 

The range of values matches those found by |Eden et al.| 
( 2015| ), varying between 0 and 80 L©/M©, with a mean value 
of 6.50 ± 0.66 L©/M© and a median value of 2.56 ± 2.06 
L©/M©. The distribution of L/M is consistent with being 
log-normal with both Anderson-Darling and Shapiro-Wilk 
tests giving probabilities of 0.107 and 0.151, respectively. 
The distribution can therefore be fitted with a Gaussian 
which has a mean of log^o(^/^) — 0-48 (i.e. L/M = 3.0) 
and a standard deviation of 1.14, giving a standard error on 
the mean of 0.075 or about 20%. These representative values 
are close to those found in an analysis of Hi-GAL data on 
infrared dark clouds in the longitude range 15° ^ ^ ^ 55° 
by Traficante et al. (2015), who found mean L/M values 
of ~3.6 and ~5.7 for starless and protostellar sources, re¬ 
spectively. The presence of a log-normal distribution tells 
us that the processes behind the distribution of L/M val¬ 
ues, both in this region and across the Galactic Plane, are 
the result of the multiplicative effect of multiple random 
processes within star-forming regions, be these the larger 
molecular-cloud environments (indicated by the log-normal 
distribution of clump-mass to cloud-mass ratios: |Eden et ah] 


2012 2013), or individual star-forming clumps. 


Eigure pT] shows the relationship between L and M for 
the subset of clumps associated with the W43 cloud. The 
gradient of the bisection-method linear regression fit to the 
distribution of points in log space is 1.76 ± 0.16. This is 
significantly steeper than the value of 1.08 ± 0.13 found by 
Urquhart et al.|p014c ), for a distance-limited set of ATLAS¬ 


GAL clumps associated with tracers of massive star forma¬ 
tion (methanol masers, Hll regions and RMS MYSOs). The 
difference may be accounted for by the inclusion of a pop¬ 
ulation of sources with low L/M values at masses below 
~10® M©, that may be relatively unevolved. 

The mean value of L/M in W43 is significantly lower 
than that of 16.0 ± 0.7 found for all ATLASGAL clumps 
associated with massive star formation but is closer to the 
6.8 ±0.7 found for methanol-maser-associated clumps alone 
by Urquhart et al. ( 2014c| . L/M is only weakly dependent 
on M, if at all, so the inclusion of higher-mass clumps in the 
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ATLASGAL sample should not account for this. The dis¬ 
tribution in Figure suggests rather that the W43 sample 
includes an excess population of sources in the mass range 
100 to 1,000 M© with low L/M values. These sources may 
be relatively young. [Nguyen- Luong et ^ (2011) found that 
the ‘future’ (i.e. incipient) SFR in W43, as traced by the 
sub-mm continuum, is around 5 times larger than the cur¬ 
rent SFR estimated from the 8-pm luminosity, and the larger 
value is consistent with the similar prediction by |Motte et| 
al. ( 2003| and see also Nguyen-Luong et al.||2013 ). 


5 SUMMARY 

Early results from the JCMT Plane Survey, part of the 
James Clerk Maxwell Telescope Legacy Survey programme 
are presented from the survey region centred on Galactic 
longitude I — 30°. The results incorporate data at about 
the 40 per cent completion stage (in both this region and 
the whole survey), the point at which the sensitivity ap¬ 
pears close to the confusion limit in this crowded area of the 
Galactic Plane. 

The contamination of the 850-pm continuum band by 
emission in the ^^CO J=3-2 line was investigated using 
COHRS data. The majority of the significant continuum 
emission has contamination levels in the region of a few per 
cent, less than the flux calibration uncertainty, and insuf¬ 
ficient to warrant the development of a detailed correction 
method at this time. 

Compact-source identification and extraction are done 
using the FellWalker algorithm resulting in a sample of 
1029 objects above a 5-cr flux threshold. The measured peak 
flux densities in this sample are compared to those in posi¬ 
tionally matched detections from the ATLASGAL 870-pm 
compact-source catalogue and found to be consistent after 
correction for the different beam sizes. 

Initial completeness tests were made using Fell- 
Walker to retrieve artificial sources from a uniform Gaus- 
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Figure 19. The relationship between embedded IR luminosity 
and clump mass for the subset of sources in the W43 JPS sample 
that have associated 70-pm detections in the Hi-GAL data. The 
red dot-dashed line indicates the linear fit to the JPS data. Also 
shown are the relationships found by ( [Urquhart et al.|2014c| > for 
ATLASGAL sources associated with: any tracer of massive star 
formation (green dashed line) and a limited-distance subsample 
(yellow dot-dash); methanol masers only (blue dotted line). The 
gradients of those lines are 1.03, 1.09 and 0.857, respectively. 



sian noise field with rms deviation equivalent to that in the 
data. The source flux densities were recovered accurately 
above peak values of 0.3 mJy per beam. Below this, fluxes 
tended to be over-estimated, perhaps due to noise-related 
biasing in the detection method. Source flux densities be¬ 
low 0.3 mJybeam”^ have been corrected for this effect. The 
completeness limit estimated in this way is of 95% complete¬ 
ness above 80 mJy beam“^, which is around 4cr in these data. 
The true completeness limit in the I — 30° region is larger 



Figure 18. The distribution of the ratio between clump (IR) 
luminosity to clump mass for the subset of sources in the W43 
JPS sample that have associated 70-pm detections in the Hi-GAL 
data. 
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than this, however, probably due to confusion. We estimate 
this level to be about 0.2 mJybeam“^. The corresponding 
limit in the integrated flux densities is ~0.8mJy which may 
be compared to a limit of 4-8 Jy in the ATLASGAL survey. 

The population of dense clumps associated with the 
W43 massive-star-forming complex clump masses is anal¬ 
ysed. The virial state of these clumps is clustered around 
the equilibrium value. The relative statistical lifetime of the 
protostellar and starless clumps (with and without a 70-pm 
detection) are about equal and the distribution of the ratio 
L/Mciump suggests that these sources are younger than the 
Milky-Way average. 
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